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Abstract 

Localization property in the disordered few-chain DNA systems with a long-range correlation is 
numerically investigated. We apply the chain system with the correlated disorder in the interchain 
and/or intrachain hoppings to the simple model of a double strand of DNA. Numerical results for 
the density of states and the Lyapunov exponent of the wave function in the two- or three-chain 
models are given. It is found that the correlation effect enhances the localization length (the inverse 
least nonnegative Lyapunov exponent) around the band center. 
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I. INTRODUCTION 



Localization phenomena in the one- dimensional disordered system have been extensively 
studied It is well known that almost all the eigenstates are exponentially localized and 
the system has a pure point energy spectrum under the presence of any disorder ^. Recently, 
the effect of the long-range correlated disorder in the potential field on the localization has 
been reported by some groups 

mm 

. The appearance of the Anderson-like metal- 
insulator transition is suggested in the strongly correlated non-stationary regime of the 
sequence jj]. As one of the realistic situation, it has been found that the base (nucleotide) 
sequence of the various genes has a strong correlation characterized by the power spectrum 
S{f) ~ (0.2 <a< 0.8) 



Moreover, the recent development o, 
of the DNA wire as a molecular device 



the nanoscale fabrication let us expect the uti 



lj,ll5[ and the realization of DNA computing 



ization 



mli7|. 



Actuall y, t he development enables us to measure the direct DNA transport phenomena 
13- Recently, Forth et al. measured the nonequilibrium current- voltage {I — V) 
characteristics in the poly(G)-poly(C) DNA molecule attached to platinum lead at room 
temperature Cuniberti et al. explained the semiconducting behavior by considering 

the base-pair stack coupled to the sugar-phosphate(SP) backbone pair 2l|. Iguchi also 
derived the semiconductivity and the band gap by using the ladder chain model of the 



double strand of DNA 



23|. In the both models, apparently, the existence of the SP 



backbone chain play an important role in the band structure due to the gap opening by the 
hybridization of the energy levels. 

On the other hand, recently Tran et al measured the conductivity along the lambda 
phage DNA (A— DNA) double helix at microwave frequencies using the lyophilized DNA in 
and also without a buffer The conductivity is strongly temperature dependent around 
room temperature with a crossover to a weakly temperature dependent conductivity at low 
temperatures. Yu and Song showed that the observed temperature dependent conductivity 
in the DNA can be consistently modeled, without invoking the additional ionic conduction 
mechanisms, by considering that electrons may use the variable range hopping for conduction 
and that electron localization is enhanced by strong thermal structural fluctuations in DNA 



2J|. Then the DNA double helix is viewed as a one- dimensional Anderson system. 



As observed in the power spectrum, the mutual information analysis and the Zipf analysis 
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of the DNA base sequence such as the human chromosome 22, the long-range structural 
correlation exists in the total sequence as well as the short-range periodicity QHQQ- 
The transport property though DNA are still controversial mainly due to the complexity of 
the experimental environment and the molecule itself. Although the theoretical explanations 
for the phenomena have been tried by some standard pict ures used in the sohd state physics 
such as polarons, solitons, electrons or holes 3, Qjisf' 



27 
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0, 31,0, 



33 



34 



35 



36|, |37[, the situation is still far from unifying the theoretical scheme. 

In the present paper we investigate the correlation effect on the localization property of 
the one-electronic states in the disordered, two-chain (ladder) [3| and three-chain models 
with a long-range structural correlation as a simple model for the electronic property in the 
DNA. The tight-binding model for the few-chain systems have the off-diagonal randomness 
as the interchain hopping integral and/or the intrachain ones. We present some numerical 
results for the density of states (DOS) and the Lyapunov exponents of the wave function. In 
particular, it is found that the correlation of the sequence enhances the localization length 
defined by using the least nonnegative Lyapunov exponent. Regardless of the parameter tun- 
ing for the numerical calculation, we would like to mainly focus on (l)suggesting the model 
and (2)giving the preliminary numerical results for the correlation effect on the localization 
in the model. 

Outline of the present paper is as follows. In the next section we introduce the models 
for DNA with a long-range correlation. In the Sect. 3 we give the results for the density of 
states. The numerical results for the Lyapunov exponent and the localization length in the 
systems are given in Sect. 4. The last section devotes for summary and discussion. 



II. MODEL 



We simplify and model the double strand of DNA by some assumptions. DNA double 
helix structure is constructed by the coupled two single strand of DNA. First, we ignore 
the twist of DNA as well as the complicated topology. (See Fig.l.) In addition to the 
simplification, we consider only the vr— electrons with the highest occupied molecular orbit 
(HOMO) states in the backbones and the base-pairs of the system. We also ignore the 
interaction between the electrons and restrict ourselves to the zero-temperature property. 
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FIG. 1: The schematic structure of a double strand of DNA. 



A. Two-chain model 

Following the basic assumptions, consider the one-electron system described by the tightly 
bin ding Hamiltonian H consisting of the two chains. This model was first given by Iguchi 



22 



23| as a model for considering the electronic properties of a double strand of DNA. 



H = Ha + Hb + Hab, (1) 
Hk = ^{kn,n\k : n){k : n\- {kn+i^n\k : n + 1) {k : n\ + h.c.)}, (2) 

n 

Hab = -Y,Vni\A:n){B -.nl + lB ■.n){A:n\), (3) 

n 

where the {\A : n >,\B : n >} denotes an orthonormalized set and (Bn+i^n) means 

the hopping integral between the nth and (n + l)th sites and An^n (-Bn,n) the on-site energy 
at site n in chain A (B), and Vn is the hopping integral from chain A{B) to chain B{A) at 
site n, respectively. The Schrodinger equation Hl"^) = becomes, 

An+l,n<Pn+l + ^n,n-10;!-l + ^„,n0^ + K.^f = E<j)^, 
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where (p^ = {A : and 



{B : Furthermore it can be rewritten in the matrix 



form, 



B 





E — Ann 




-An — lri 





-^n + ln 


Bn + ln 






Vn 


E — Bnn 









Bn+ln 


Bn+ln 






1 











v 





1 









/ 6^ \ 



J 



Tin) 



(4) 



Generally speaking, we would like to investigate the asymptotic behavior {n oo) of the 
products of the matrices M„ = n^T(/c). According to the parameter sets given by Iguchi 
we set An+i^n = -Bn+i,n = a(= at odd (even) site ra, respectively and K = at even 
sites (phosphate sites) for simplicity. The chain A and B are constructed by the repetition 
of the sugar-phosphate sites, and the inter-chain hopping K„ at the sugar sites come from 
the nucleotide base-pairs, i.e., A — ToiG — C pairs. (See Fig.2(a).) 




FIG. 2: Models of the double starnd of DNA. (a) the two-chain model, (b) the three-chain model. 



B. Three-chain model 

When we pay attention to the overlap integral between the nucleotide base-pairs, the 
two-chain model can be easily extended to the three-chain one described by following Hamil- 
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tonian. (See Fig.2(b)0 



H = Ha + Hb + He + Hac + Hbc-, 
Hac + Hbc = -Y,^Vn\A:n){C ■.n\ + Un\B:n){C ■.n\ + h.c.] 



(5) 
(6) 



where the {\A : n >^\C : n >^\B : n >} denotes an orthonormahzed set and C^+i^n means 
the hopping integral between the nth and {n + l)th sites and Cn,n the on-site energy at site 
n in chain C, and Vn and f/„ are the hopping integral between the chains. The Schrodinger 
equation H\^) = becomes, 



Cn+l,n0^+l 



a 



ill + An,n<l)n + Vn<Pn = E(j)t 



A 



7i,n-10n-l + Cn,n4'n + ^nS^n 



U, 



C 



Bn,n-l(t>n-l 



n 



where 



kC 



Bn,n<P. 

{C : n|\l/). It can be rewritten in the matrix form 



^ n. 



c 

n'Tn 



T n) 



^n— 1 



,T{n) 




(7) 



where $„ = 0„ , 0„ ) and / and O denote the three-dimensional unit and zero matrices, 
respectively. The matrices Hd and Ho are given as follows: 

/ E-A 

nn ^n— In 



V 



E — Ann 


Vn 







E — Cnn 


C/n 




C'n + ln 


C*n+ln 







E — Bn„ 


Bn + ln 


Bn+ln 



, Ho 









Cn — ln 



\ 
















En— In 
En-\-ln / 



(8) 



(See 



In this paper we set A„+i „ = -B„+i,„ = a(= 6) at odd(even) site n, and Vn = Un- 
Fig. 2(b).) In addition, in the numerical calculation we set the on-site energy as Ann = Bnn = 
Cnn = for simplicity. The sequence {C„„+i} can be also generated by corresponding to 
the base-pairs sequence. The localization properties in the simple few-chain models with 
the on-site disorder have been extensively investigated jsOj. 



C. The correlated sequence 



The correlated binary sequence {V^} and /or jCnn.-t-i I of the hopping integrals can be 
generated by the modified Bernoulli map P, 3, 
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6 



( 



(9) 



where Jq = [0, 1/2), Ji = [1/2, 1). 5 is a bifurcation parameter which controls the correlation 
of the sequence. 

In the ladder model we use the symbolized sequence {Ki} by the following rule as the 
interchain hopping integral at odd sites n: 



In the numerical calculation, we use Wqc = ^at/'^ for simplicity. Then the artificial binary 
sequence can be roughly regarded as the base-pair sequence as observed in the A— DNA or 
the human chromosome 22. The correlation function C{n){=< VnoVn^+n >) {no = 1, n is 
even ) decays by the inverse power-law depending on the value B as C{n) ~ n~^-^ for large 
n (3/2 < B < 2). The power spectrum becomes S{f) ~ f'^^ for small /. We focus 
on the Gaussian (1 < B < 3/2) and the non-Gaussian stationary region 2/3 < 5 < 2, 
corresponding to the DNA base-pairs sequence (0.2 < a < 0.8). 

In the three-chain model we assume that the phosphate sites are renormalized in the 
hopping integral between the sugar sites for simplicity. The interchain hopping integrals = 
Un can be generated by the same way to the two-chain model for every site n. Furthermore, 
we use the successive sequence {Xn, Xn+i} when we make the correlated binary sequence 
{Cnn+i} as the hopping integral of the middle (nucleotide) chain as follows: 




Wat^WtA (Xnelo) 

Wgc = WcG {Xn e h). 



(10) 




nn+1 ^ 



WaT-GC {Xn e lQ,Xn+l G h) 



(11) 



Wgc -AT {Xn e h,Xn+i e h) 
^ Wgc-gc {Xn e h,Xn+i eh). 
In the numerical calculation we assume the following rules for simplicity. 



Wgc-gc — Wat -at/ 2 
< Wat-gc — (Wat- AT + Wgc-gc)/ 2 
Wgc -AT — Wat-gc- 
As a result, the parameters are Wat and Wat -at- 



(12) 
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III. DENSITY OF STATES AND LYAPUNOV EXPONENTS 



In this section, we show the DOS and the Lyapunov exponents for the periodic and 
disordered two-chain model in order to give some notations and to confirm the numerical 
reliability. 

Let us consider the periodic case to confirm the gap opening mechanism due to hy- 
bridization by the interchain coupling. Figure 3 shows the DOS for various periodic cases. 
As seen in the Fig. 3(b) the energy spectrum consists of four bands when Wat = Wgc and 
Ann = Bnn = 0. The band gap Eg{v) at the center is given by Eg{v) = + f^) — v, 

where Eg = 2{a — b) means the band gap for the single strand of DNA. The other band gap 
Ag{v) = V + ^[y/{Eg + 1>^) + ■\/(4(a + by + v"^)] appears in between the lowest and the lower 
middle bands where v > 2abl ^{a^ + 6^). Iguchi showed the semiconductivity of the double 
strand of DNA based on the band structure 0, 13- The energy spectrum in the case a = b 
is given in Fig.3(a). The band gap at the center disappears. Furthermore, the energy band 
for case Wat 7^ WcciWcc = Wat^^) is shown in Fig.3(c). Each band splits and the total 
spectrum consists of eight bands. 

Figure 4 shows the DOS as a function of the energy for the binary disordered system. 
The sequence of the interchain hopping takes an alternative value Wqc ot Wat- It is found 
that some gaps observed in the periodic cases close due to the disorder corresponding to the 
base-pairs sequence. 

It is noted that in the bipartite lattice with an even number of sites the energy spectrum 



is striking symmetric around E = 



43 



44 



45 



46 



3- 



47,48^. In this case, related to the chiral 



universality class, the E = states are non-localized in any dimension. However, the states 
with the Dyson singularity 3] at the center disappear when the randomness is introduced 
in the on-site energy at the backbone chains. We comment the point in the next section 
again. 

We give a preliminary numerical result of the energy dependence of the Lyapunov expo- 
nents. The definition is given by, 

7i= lim log a, (MtM„), (13) 

where crj(...) denotes the ith eigenvalue j41j. As the transfer matrix T{n) is symplec- 
tic, the eigenvalues of the M\Mn have the reciprocal symmetry around the unity as 
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FIG. 3: DOS as a function of energy for the periodic cases in the two-chain model. {a)WAT 
Wgc = 2.0, a = 1.0,6 = 1.0. {h)WAT = Wgc = 2.0, a = 1.0,6 = 0.5. {c)Wat = 2.0, Wgc 
1.0, a = 1.0, h = 1.0. The on-site energy is set at Ann = Bnn = 0. 



=71 



old 



-Id 



-71 



where 71 > 72 > • • • 7^ > 0. The d denotes the number of channels, 



i.e. d = 2 in the two-chain model and d = 3 in the three-chain model. 

Figure 5 shows all of the Lyapunov exponents 7j(z = 1, 2, 3, 4) as a function of the energy 
to confirm the reciprocal symmetry. The energy regions where two exponents are positive 
(71 > 72 > 0) correspond to the energy gap in Fig. 3 except for the vicinity of the singular 
point E = {). The zero Lyapunov exponent corresponds to the brunch cut in the analyticity 
of the energy spectrum. 

The Lyapunov exponents are related to the DOS p{E) as an analogue of the Thouless 
relation as 



apunov ( 



d » 

J^7i(E)~ \n\E~E'\p{E')dE'. 



(14) 



Note that the Thouless relation {d = 1) is true only for the positive Lyapunov exponents. 
One can see that the singularity of the largest Lyapunov exponent is strongly related to the 
singularity of the DOS. 

Furthermore, it is found that for the thermodynamic limit the largest channel-dependent 
localization length C,d = l/7d determines the exponential decay of the Landauer conductance 
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C/3 

o 

a 




FIG. 4: DOS as a function of energy for the binary disordered cases in the two-chain model. 
{a)WAT = 2.0, Wgc = 1-0, a = 1.0,6 = 1.0. (h)WAT = 2.0,Wgc = 1-0, a = 1.0,6 = 0.5. The 
on-site energy is set at Ann = Bnn = 0. 




FIG. 5: Lyapunov exponents as a function of energy for the periodic cases in the two-chain model. 
{a)WAT = Wgc = 2.0, a = 1.0,6 = 1.0. {h)WAT = 2.0, Wgc = 1-0, a = 1.0,6 = 1.0. The on-site 
energy is set at Ann = Bnn = 0. 

g which is in units of /h at zero temperature and serves as the locahzation length of the 
total system of the coupled chains ^3, ^ ■ 

d 



. , cosh^ 



2m 

exp --— — . 



(15) 



Recently, electron transport for the molecular wire between two metalic electrorodes has 



been also investigated by several techniques j5J, |55 , 



56|,|57| 
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IV. CORRELATION EFFECTS 



In this section, we consider the correlation effect on the locahzation property of the 
disordered case by using the modified Bernoulh model. Figure 6 shows some eigenstates 
for the correlated case, B = 1.8. Apparently the eigenstates are exponentially localized. In 
particular, the energy dependence of the Lyapunov exponents and the localization length 
are investigated in order to characterize the localization. We used a large sample with the 
system size = 10^ for the numerical calculation. 




site n 



FIG. 6: Some eigenfunctions for the binary correlated disordered cases B = 1.8 in the two-chain 
model (A^ = 256). The amplitude of the wave function for A and B chains are over- written. 
Wat = 2.0, Wgc = 1-0, a = 1.0, b = 0.5. The on-site energy is set at Ann = Bnn = 0. 

A. Two-chain model 

First, we show the expectable results by a simple analytical calculation. We can write 
the Schrodinger equation in the following matrix form when the hopping integral and the 
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on-site energy of the backbone chains are the same values (a = b, Ann = Bnn), 



^tl +^tl\ 1 ( E- Ann -Vn 

+ €-1 



-Vn E — Br, 





1 















It is easily seen that the equations can be decoupled to the one-dimensional Anderson model 
with diagonal randomness by the unitary transform U as follows; 



(16) 





where 



V'f+l+V'f-ij M '^n E + Vn 




(17) 



f^n — '^{Ann + Bnn)i — '^{Ann " Bnn)- (l^) 




-2-10 1 

E 



FIG. 7: Lyapunov exponents 7i(i = 1, 2) as a function of energy in the correlated two-chain model. 
The parameters Wat = 2.0, a = 1.0, b = 1.0. The on-site energy is set at Ann = Bnn = 0. The 
number of the matrix product is = 10®. 

In our case, /i„ = Ann, '^n = 0. Accordingly, the behavior of the Lyapunov exponents 
are essentially the same to that in the one-dimensional system with the sublattice diagonal 
disorder. After all the Vn is non-zero only if n is odd, i.e. V2„-i is random and V2„ = 0. 
Sometime the model is called the "periodic random alloy" or the "dilute Anderson model" 



43L I46| . Although in our case the energy at the deterministic even sites are zero, the case 

12 



0.20 ■ 



0.15 ■ 



0.10 ■ 



0.05 ■ 



0.00 - 




FIG. 8: Lyapunov exponents 7i(i = 1, 2) as a function of energy in the correlated two-chain model. 
The parameters Wat = 2.0, a = 1.0, b = 0.5. The on-site energy is set at Ann = Bnn = 0. The 
number of the matrix product is = 10^. 



where the deterministic sites are non-zero but constant is trivially obtained by shifting the 
energy. The extended state at = appears in the one-dimensional periodic random binary 



alloys due to a resonance |46|, l58l |. It is demonstrated that the occurrence of the extended 



states in the vicinity of the site energy of the deterministic part, i.e. E = in our case. The 
delocalization can be understood in terms of the long-range correlation due to the perfect 
periodicity in the deterministic sites. 

In addition, we introduce another long-range correlation due to the base-pair sequence on 
the odd sites V2n-\- Figures 7 and 8 show the energy dependence of the Lyapunov exponents 
(71 and 72) in some correlated cases [B = 1.0, 1.5, 1.8). It found that the correlation enhances 
the localization length 7^^) around \E\ < 1, although the largest Lyapunov exponent 
7i is almost remained. 

In Fig. 9(a) the localization length = ^/id defined by the least nonnegative Lyapunov 
exponent is shown. The electronic states whose energy is close to the resonance turn out 
to be extended in the sense that the localization length is the same order or larger than 
the system size. Accordingly the correlation effect is not so strong on the global feature in 
the energy dependence of the density of states, particularlly the correlation enhances the 
localization length in the energy regime \E\ < 1. 
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FIG. 9: Localization length I/72) as a function of energy in the correlated (a)two-chain model 
and (b)three-chain model. The parameters Wat = "^-OtWat-AT = 1-0, a = 1.0,6 = 0.5. The 
on-site energy is set at Ann = Bnn = 0. The number of the matrix product is = 10^. 

B. Three-chain model 

In this subsection, we confirm the numerical results for the three-chain model when 
compared to that in the two-chain model. It must be noted that if we take into account of 
a vr — vr interaction between the stacked base pairs, the model system becomes three chains 
model. (See Fig.2(b).) 

As has been done in the two-chain case we can rewrite the Schrodinger equation in the 
matrix form as follows. 









\ 

1 


(e 


_ A 


-Vn 





\ 


( 




f 


1>^^ 




( 






1 




V„ 

tc 


E — Cnn 
tc 


tc 






= R 


( 




(19) 


V 


+ 






\ 





Un 


E - B 


I 






V 


) 





We set a = 6 = 1, C„,n+i = tc for simplicity. It should be noted that the matrix R can 
not be diagonalized by any unitary transform because the matrix is not a normal matrix, 
R^R 7^ when tc 7^ 1 (Toeplitz theorem). Accordingly the simple expectation in the 

previous subsection is not suitable for the three-chain model. 

Figures 10 and 11 show the energy dependence of the Lyapunov exponents in the three 
chain cases. We can observe that all the Lyapunov exponents 7j(i < d) are changed by the 
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correlation. The least nonnegative Lyapunov exponent 73 is diminished by the correlation. 
The localization length are shown in Fig.9(b). We can see that the localization lenght 
diverges at the band center E = 0. 

Generally speaking, in the quasi-one-dimensional chain with the hopping disorder the 
singularity of DOS, the localization length and the conuctance at the band center depend 



on the parity, bipartiteness and the boundary condition |46l . |48| . As the discussion about 
the details is out of scope of this paper, we give the simple comment. Note that the parity 
effects appear in the odd number chain with the hopping randomness. In the odd number 
chain with a hopping randomness, only one mode at = is remained as the extended 
state, i.e. 7^ = 0, and the other exponents are positive, ■jd-i > ■ ■ ■ > 71 > 0. The behavior is 
seen in Fig. 10 and 11. Then, the non-localized states with 7 = determine the conductance. 
Although we ignored the bipartite structure in the three-chain model for simplicity, if we 
introduct the bipartiteness in the intrachain hopping integral Vn{= Un) the other delocalized 
state due to the chiral symmetry appears at E' = as seen in the ladder model. 

Here, we stress again that in the more realistic model of DNA the parity effects and the 
chiral symmetry are broken down by the on-site energy fluctuation. Indeed, Fig. 12 (a) shows 
the energy dependence of the Lyapunov exponents in the system which the on-site energy 
Arm{= Bnn) in the side chains is also randomly genarated as well as the hopping disorder 
in Fig. 11. In this case we can confirm that the least nonnegative Lyapunov exponent, 73, is 
positive. It is found that the extended state at = disappears and the localization length 
becomes finite. (See Fig. 12(b).) As a result, the localization properties can be changed by 
the introduction of the correlation and the fluctuation of the sequence even in scope of our 
models. 



V. SUMMARY AND DISCUSSION 

We numerically investigated the correlation effect on the localization property of the 
one-electronic states in the disordered two-chain (ladder) and three-chain models with the 
long-range structural correlation as the simple models for electronic property in the double 
strand of DNA. The results are summarized as follows. 

(1) We gave the tight-binding models for the few-chain systems with the off-diagonal 
randomness as the interchain and/or the intrachain hopping integrals. 
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FIG. 10: Lyapunov exponents 7j(i = 1,2,3) as a function of energy in the correlated three-chain 
model. The parameters Wat = 2.0, Wat-AT = 1-0, a = 1.0,6 = 1.0. The on-site energy is set at 



(2) A simple interpretation for the results obtained in the two-chain model was given by 
means of the unitary transform. 

(3) The correlation enhances the localization length (7^^) around \E\ < 1, although the 
71 is almost remained, in the two-chain model. 

(4) In the three-chain model, all the Lyapunov exponents are changed by the correla- 
tion effect, particularlly, the least nonnegative exponent 73 is diminished by the long-range 
correlation. 

(5) The divergence of the localization length around the E = disappears by the intro- 
duction of the fluctuation of the on-site energy. 

The details of the energy dependence, the correlation dependence and the potential 
strength dependence on the DOS and the Lyapunov exponents, will be given elsewhere 



Although we concentrated on the localization property of the correlated system for the 
double strand of DNA by using simple values for the parameters, the more realistic Hamilto- 
nian matrix elements are obtained from the oligomer calculation using the extended Hiickel 



In the experiment of the conductance property of the DNA both the temperature de- 








theory with the dementi's parametization 
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FIG. 11: Lyapunov exponents 7j(i = 1,2,3) as a function of energy in the correlated three-chain 
model. The parameters Wat = 2.0, Wat-AT = 1-0, a = 1.0,6 = 0.5. The on-site energy is set at 

-^nn — — Cnn — 0. 

pendence and the temperature effect are important. Finite temperature can also reduce the 
effective system size and leads to the changes in the transport property. Moreover, the effects 
of the stacking energy and of temperature can be considered by introducing the fluctuation 
in the hopping energy such as the Su-Schriefer-Heegar model for polyacetilene [g^. It is 
interesting to investigate the quantum diffusion in order to reveal the conducting properties 



of electrons, polarons and solitons in the fluctuating ladder models [59[. In the present 
model, although we used the binary correlated sequence, the four letter virsion (A,T,G,C) 
is also interesting An extension for the four letter virsion with a long-range correlated 
hopping disorder has been constructed 
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FIG. 12: (a) Lyapunov exponents ^{{i = 1,2,3) as a function of energy in the disordered three- 
chain model (B = 1.0). (b) The localization length defined by = 1/73- The parameters 
Wat = 2.0, W^T-AT = 1.0, a = 1.0, b = 0.5 and C„„ = 0. The on-site energy of the side chains, 
Ann{= Bnn), IS randomly generated. 

study. 
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